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including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

APR  2013 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2013  to  00-00-2013 

4.  TITLE  AND  SUBTITLE 

Maximum  Likelihood  Detection  of  Low  Rate  Repeat  Codes  in  Frequency 
Hopped  Systems 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

MITRE  Corporation,1155  Academy  Park  Loop,CoIorado 
Springs,CO,80910-3761 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

In  time  diversity  wireless  and  satellite  communication  multiple  copies  of  the  same  signal  segment  are 
transmitted  during  different  time  intervals  to  improve  signal  detection.  If  the  communication  system  is 
frequency  hopped,  coherent  detection  is  infeasible.  In  the  traditional  approach  to  deal  with  this  problem, 
the  receiver  uses  only  the  reference  signals  for  phase  shift  correction  and  combines  soft  symbol  decisions 
obtained  independently  from  each  copy.  In  this  paper,  we  develop  the  corresponding  theoretical  maximum 
likelihood  (ML)  detection  problem,  its  solution,  and  a  computationally  efficient  algorithm  that  is  an 
approximate  ML  solution.  We  present  several  simulation  experiments  and  results.  The  experiments  include 
phase  drifts  allowable  in  practical  systems  and  randomized  variations  in  the  locations  of  reference  symbols. 
Results  indicate  power  savings  of  up  to  2.0  decibel  (dB)  over  the  traditional  method  for  different  system 
configurations.  They  also  show  that  for  short  data  segments  used  in  systems  over  severely  degraded 
channels,  the  results  from  the  theoretical  solution  and  our  algorithm  are  virtually  indistinguishable.  The 
main  impact  of  this  research  is  that  survivable  and  protected  communication  systems  can  take  advantage 
of  our  new  signal  combining  algorithm  that  offers  considerable  power  savings. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

1 9a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

9 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


This  page  intentionally  left  blank. 


Approved  for  Public  Release;  Distribution  Unlimited.  13-1340 

Maximum  Likelihood  Detection  of  Low  Rate  Repeat 
Codes  in  Frequency  Hopped  Systems 


G.  R.  (Datta)  Dattatreya 
The  MITRE  Corporation 
Colorado  Springs 


Abstract — In  time  diversity  wireless  and  satellite  communication, 
multiple  copies  of  the  same  signal  segment  are  transmitted  during 
different  time  intervals  to  improve  signal  detection.  If  the 
communication  system  is  frequency  hopped,  coherent  detection  is 
infeasihle.  In  the  traditional  approach  to  deal  with  this  problem,  the 
receiver  uses  only  the  reference  signals  for  phase  shift  correction 
and  comhines  soft  symbol  decisions  obtained  independently  from 
each  copy.  In  this  paper,  we  develop  the  corresponding  theoretical 
maximum  likelihood  (ML)  detection  problem,  its  solution,  and  a 
computationally  efficient  algorithm  that  is  an  approximate  ML 
solution.  We  present  several  simulation  experiments  and  results. 
The  experiments  include  phase  drifts  allowable  in  practical  systems 
and  randomized  variations  in  the  locations  of  reference  symbols. 
Results  indicate  power  savings  of  up  to  2.0  decibel  (dB)  over  the 
traditional  method  for  different  system  configurations.  They  also 
show  that  for  short  data  segments  used  in  systems  over  severely 
degraded  channels,  the  results  from  the  theoretical  solution  and  our 
algorithm  are  virtually  indistinguishable.  The  main  impact  of  this 
research  is  that  survivable  and  protected  communication  systems 
can  take  advantage  of  our  new  signal  combining  algorithm  that 
offers  considerable  power  savings. 

Keywords — diversity  communication;  joint  phase  estimation  and 
detection;  maximum  likelihood;  signal  combining;  frequency 
hopping;  repeat  codes 

I.  INTRODUCTION 

Satellite,  airborne,  and  terrestrial  wireless  communications 
through  severely  degraded  and/or  contested  channels  require  the 
use  of  low  rate  coding.  A  practical  approach  to  achieving  such 
low  rate  codes  is  to  first  construct  forward  error  correction 
(EEC)  codes  of  half  or  one -third  rate  and  then  transmit  multiple 
copies  of  these  code  words.  Protected  communication  systems 
employ  frequency  hopping  and  the  dwell  time  of  each  carrier 
frequency  is  too  short  for  carrier  phase  acquisition.  This  prevents 
coherent  demodulation  [1]  at  the  receiver.  The  symbol  sequence 
transmitted  over  a  dwell  time  is  called  a  hop.  A  hop  includes 
reference  symbols  of  known  values  in  addition  to  the  data 
symbols.  In  severely  degraded  channel  conditions,  the  system 
uses  larger  and  larger  symbol  transmission  times  to  increase  the 
signal  to  noise  ratio  (SNR).  Under  such  conditions,  a  hop  would 
contain  only  a  modest  number  (for  example,  a  few  tens)  of  data 
symbols  and  a  few  reference  symbols.  Transmitted  hops  are 


subjected  to  unknown  phase  shifts  during  radio  propagation. 
Received  hops  should  be  demodulated  and  their  soft  decisions 
input  to  the  decoder.  Typically,  the  demodulator  uses  a  local 
oscillator  with  an  arbitrary  phase  and  outputs  a  complex  valued 
number  for  each  symbol.  We  refer  to  each  of  these  complex 
numbers  (corresponding  to  a  symbol)  as  a  signal.  The  rest  of  the 
receiver  should  be  designed  to  deal  with  the  unknown  arbitrary 
phase  shift  in  the  sequence  of  such  signals  produced  by  the 
demodulator.  Noise  on  reference  symbols  introduces  errors  in 
the  estimation  of  such  phase  shift.  It  is  important  to  minimize 
this  effect  through  the  proper  use  of  all  of  the  received  reference 
and  data  signals  in  a  hop. 

Multiple  copies  of  a  hop  are  subjected  to  independent  unknown 
phase  shifts  and  additive  noise.  At  a  first  glance,  it  might  appear 
that  the  best  approach  for  the  receiver  to  deal  with  two  copies  is 
to  apply  the  ML  phase  shift  and  detection  estimation  algorithm 
due  to  Mackenthun  [2]  separately  for  each  of  the  two  hops, 
obtain  the  sequences  of  symbol  or  bit  log  likelihood  ratios 
(LLRs),  and  then  add  the  LLRs.  The  following  argument  shows 
that  such  an  approach  is  suboptimal.  When  we  apply  the  ML 
algorithm  separately  on  the  two  copies,  the  algorithm  allows  for 
the  detected  data  symbols  at  a  particular  position  (say  the 
seventh  position  in  the  two  hops)  to  be  different.  Clearly,  this  is 
suboptimal  and  a  good  algorithm  should  constrain  that  the 
symbol  detections  at  the  corresponding  positions  of  the  two  hops 
be  the  same.  In  this  paper,  we  study  joint  ML  phase  estimation 
and  detection  of  multiple  copies  of  a  hop. 

Some  background  and  related  literature  are  reviewed  in  Section 
II.  Section  III  is  devoted  to  the  development  of  the  ML  signal 
combining  problem  and  its  solution.  The  global  objective 
function  for  c  copies  of  an  M-ary  Phase  Shift  Keyed  (MPSK) 
hop  with  n  data  symbols  and  r  reference  symbols  in  each  copy  is 
developed  and  simplified.  Complexities  of  the  theoretical  ML 
solution  are  discussed.  Our  signal  combining  algorithm,  which  is 
an  approximation  to  the  ML  solution,  is  developed  in  Section 
IV.  Section  V  presents  several  simulation  experiments  and 
results.  The  experiments  include  different  burst  modes  (different 
numbers  of  data  and  reference  symbols  in  a  hop),  different  types 
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of  modulation,  phase  drift  caused  by  local  oscillator  frequency 
error,  and  randomized  locations  of  reference  symbols  within  a 
hop.  Section  VI  concludes  the  paper  with  highlights  of  the 
results  and  the  impact. 

II.  BACKGROUND  AND  RELATED  LITERATURE 
Consider  a  coherent  communication  system  [1]  that  receives  two 
noisy  signal  segments  corresponding  to  identically  transmitted 
copies.  Let  Eg  be  the  energy  per  data  symbol  at  the  receiver  in 
each  copy  and  let  N-^  and  N2  be  the  spectral  densities  of  the 
additive  white  Gaussian  noise  (AWGN)  corrupting  the  two 
copies,  respectively.  If  the  corresponding  noise  standard 
deviations  are  Ui  and  cr2,  combining  the  two  signal  segments  is 
through  weighted  averaging,  with  weights  proportional  to  —  and 

—  (Chase  [3]).  If  =  (T2,  the  process  of  combining  is  simply 
02 

averaging  the  two  signal  sequences  before  obtaining  soft 
decisions.  In  this  case,  the  effect  of  transmitting  two  copies  of  a 
data  segment  is  the  equivalent  to  having  transmitted  only  the 
original  segment  (and  no  copy)  at  double  (or  3  dB  above)  the 
power.  As  mentioned  earlier,  wireless  and  satellite 
communication  systems  that  use  frequency  hopping  cannot 
implement  coherent  demodulation,  and  the  transmitted  hops 
include  known  reference  symbols  to  help  with  demodulation  and 
to  obtain  soft  decision.  Consider  MPSK  systems  with  n  data 
symbols  and  r  reference  symbols  in  each  hop.  Assume  that  all 
the  reference  symbols  are  transmitted  with  a  zero  phase.'  The 
simplest  approach  to  estimate  the  phase  shift  is  to  average  the  r 
received  reference  signals.  A  simple  way  to  imagine  the  joint 
maximum  likelihood  phase  estimation  and  data  symbol  detection 
is  to  construct  the  likelihood  function  for  each  of  all  the  M” 
possible  symbol  sequence  detections  and  select  the  one  with  the 
highest  likelihood  function.  If  necessary,  we  can  then  obtain  the 
best  estimate  of  the  phase  shift  by  assuming  that  the  best 
sequence  detection  is  correct.  An  alternative  approach  for  joint 
ML  detection  and  estimation  is  to  consider  all  possible  values  of 
phase  shifts,  evaluate  the  sequence-detection  and  the 
corresponding  likelihood  function  for  each  case  of  the  phase 
shift,  and  then  select  the  best  sequence  detection.  In  this 
procedure,  as  we  increase  the  candidate  phase  shift  variable, 
corresponding  symbol  decision  boundary  radii  rotate.  When  (and 
only  when)  a  boundary  crosses  a  data  signal,  the  detection  of 
only  that  symbol  changes.  This  induces  a  change  in  the  behavior 
of  the  likelihood  function  at  every  value  of  the  phase  shift 
variable  at  which  a  data  signal  falls  on  one  of  the  decision 
boundary  radii.  The  continuous  variation  of  phase  shift  variable 
between  two  successive  decision  boundary  crossings  is 
irrelevant.  This  leads  to  an  efficient  algorithm  of  computational 
complexity  nlogn  developed  by  Mackenthun  [2]  and 
independently  by  Sweldens  [4].  Motedayen-Aval  and 
Anastasopoulos  [5]  generalize  this  result  and  identify  a  class  of 
problems  for  which  the  hard  or  soft  symbol  sequence  decisions 


'in  practice,  each  reference  symbol  can  be  any  valid  symbol.  The  receiver  is 
required  to  know  the  exact  value  of  the  reference  symbol  used  at  each  location. 
The  choice  of  the  reference  symbol  does  not  affect  the  BER. 


in  the  presence  of  unknown  parameters  can  be  exactly  evaluated 
with  only  polynomial  complexity  with  respect  to  the  sequence 
length. 


III.  ML  SIGNAL  COMBINING  WITH  UNKNOWN 
PHASE  SHIFTS 

Consider  the  problem  of  sequence  detection  from  two  copies  of 
a  hop  received  with  unknown  and  independent  phase  shifts  and 
AWGN.  The  joint  probability  density  function  we  need  to 
maximize  for  signal  combining  of  two  copies  is 


^2(n+r)^5n+r)  2(n+r) 


exp 


i  =  l 

jdi  _  pjXdi+(p2)\ 


(1) 


In  the  above  probability  density  function,  the  sequences 
and  {hi, , hn+r)  the  magnitude  and  phase 
parts  of  the  first  copy  of  the  received  hop.  The  sequences 
(c^,  , Cn+r)  {di,  , dn+r)  ^Tc  those  of  the  second  copy  of 

the  hop.  cr^  and  cr|  are  the  variances  of  the  complex  AWGN  in 
the  two  hops,  respectively.  is  the  phase  angle  of  the  data 
symbol  i.  If  these  two  variances  are  different,  we  can  normalize 
the  magnitudes  of  the  signals  in  the  received  hops.  Therefore,  we 
can  assume  without  loss  of  generality  that  cr^  =  cr|  =  cr^.  The 
phase  angles  cpi  and  (p2  are  the  unknown  phase  shifts  in  the  two 
hops.  As  earlier,  there  are  r  reference  symbols  in  each  hop.  The 
optimization  problem  becomes 


min 

<Pi,<P2,{e} 


+ 


_  gjCSi+cpt) 


n+r 


i=l 


—  e/C0i  +  <P2) 


(2) 


where  0  is  the  sequence  61, ...  of  phase  angles  for  the  n  data 
symbol  decisions  and  {6}  is  the  set  of  phase  angle  sequences  of 
all  possible  data  sequence  decisions.  With  this  background  on 
two  copies  of  a  hop,  we  can  generalize  to  multiple  copies.  Now, 
let  there  be  c  copies  of  hops.  Let  each  hop  contain  n  data 
symbols  and  r  reference  symbols  for  a  total  of  n  -f  r  symbols 
per  hop.  The  order  of  interleaving  the  data  symbols  and 
reference  symbols  within  the  hop  is  irrelevant.  We  just  need  to 
be  able  to  rearrange  them  at  the  receiver  to  restore  the  same 
order  for  both  the  hops.  Therefore,  we  assume  without  loss  of 
generality  that  the  sequences  of  transmitted  symbols  in  all  the 
copies  of  a  hop  are  identical.  Let  the  magnitude  of  the  i-th  data 
signal  in  the  y'-th  received  copy  of  the  hop  be  a/y  and  the  phase 
angle  of  the  same  be  hjy.  The  unknown  phase  shift  of  the  y-th 
copy  is  (pj.  The  magnitude  of  the  i-th  reference  signal  in  the  y'-th 
copy  of  the  hop  is  Sjy  and  the  phase  angle  of  the  same  is  tjy.  We 
can  express  the  joint  density  function  as  earlier  and 
mathematically  manipulate  the  expression  that  is  required  to  be 
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maximized.  After  eonsiderable  manipulation,  we  obtain  the 
following  expression  for  optimization. 


max /(a,  b,  s,t,(p,  9) 


j=l  Lvi=i 

ir  > 

^S^jCOs{(pj  -  ty)j 


(3) 


On  the  left  hand  side  (LHS)  of  the  above  equation,  the  bold 
letters  represent  the  sequenees  (or  vectors)  corresponding  to  the 
subscripted  variables  found  on  the  right  hand  side  (RHS).  This 
optimization  problem  has  c  continuous  phase  angle  variables, 
(Pj,  each  varying  over  [0,  2n)  radians,  and  n  discrete  variables 
6i,  each  of  which  takes  one  of  the  M  phase  angles  in  the  MPSK 
system. 

As  in  the  case  of  optimization  for  only  one  copy  of  a  hop,  for 
any  candidate  phase  vector  parameter  <p,  there  exists  a 
correspondingly  best  sequence  of  symbols  decisions,  6*  {<p)  and 
these  are  straightforward  to  evaluate,  separately,  for  each  data 
signal  at  index  i  =  1,  ,  n.  That  is,  8*  is  obtained  by  optimizing 

as  follows: 

C 

max  /j {(p,  Bi)  =  m^ax ^  cos(  cpj  -  b^j  +  .  (4) 

;=i 


Note  that  if  the  phase  shift  parameter  (p  is  given,  the  reference 
signals  do  not  influence  the  symbol  detections.  Alternatively,  if 
we  are  given  a  sequence  of  symbol  decisions  0,  we  can  evaluate 
the  correspondingly  best  reference  phase  vector  (p*,  with 
components  (pj,  j  =  1,  ,  c.  The  objective  function  with  known 

9i  values  and  unknown  (pj  is  mathematically  manipulated  and 
this  results  in 


f((p)  =  ^  |cos((py)  ^  aijCos{bij  -  ^i)  +  ^  N;  cos(tiy) 


Let 


;=i 


+  sin((py) 


IL 

^  aijsm{bij  -  Bi) 


+ 


^  Sij  sin(tjy) 


^;  = 

y  ajyCOs(fyy  -  -f  y  Sij  cos(tjy) 

.  (6) 

i=l  i=l 

■  n  r 

Bj  = 

^  aysin(fyy  -  ^i)  +  ^  N;  sin(fy) 

■  i=l  i=l 

,  (7) 

Aj  +  Bf, 

(8) 

(5) 


and  j/^yuniquely  defined  by 


B 

sin(i/)y)  = -r,and  (9) 

W 

cos(i/)y)  =  (10) 

We  now  have 

C 

fi<p)  =  ^  Cj  cos{(pj  -  ipj).  (11) 

/=! 

Maximizing  fi<p)  is  now  very  simple.  The  solution  is  given  by 

<Pj=^j>  j  =  l,-,c.  (12) 

Now,  consider  the  behavior  of  the  objective  function  over  the 
range  of  all  possible  c  phase  shift  variables,  the  c-dimensional 
hypercube  with  each  edge  of  the  cube  extending  over  [0,27r). 
This  hypercube  is  implicitly  partitioned  into  many  sub-regions 
such  that  the  best  sequence  detection  is  the  same  for  all  phase 
angle  vectors  in  a  given  sub-region.  The  boundaries  of  these  sub- 
regions  are  nonlinear  (as  opposed  to  being  hyper-planes).  The 
objective  function  is  continuous  everywhere  in  the  hypercube, 
but  it  is  not  differentiable  at  points  on  the  hyper-surface 
boundaries  separating  these  sub-regions.  In  general,  some  of 
these  (at  least  one)  sub-regions  have  a  point  each  at  which  all 
partial  derivatives  of  the  objective  function  vanish.  If  we  start 
with  any  point  in  such  a  sub-region,  evaluate  the 
correspondingly  best  data  symbol  detections,  and  re-estimate  the 
phase  shift  corresponding  to  these  detections,  the  result  is  the 
point  with  zero  partial  derivatives  and  a  local  optimum.  If  we 
start  from  any  point  in  a  sub-region  that  does  not  have  a  point 
with  zero  partial  derivatives,  evaluate  the  correspondingly  best 
data  symbol  detections,  and  re-estimate  the  phase  shift  vector 
based  on  the  detections,  we  would  have  moved  out  of  the  sub- 
region  into  a  different  sub-region.  This  can  be  confirmed  by  re¬ 
evaluating  the  best  detections  for  the  re-estimated  phase  shift 
vector  and  comparing  it  with  the  detections  in  the  first  iteration. 
Proofs  of  these  properties  of  the  objective  function  follow  as 
extensions  to  corresponding  properties  in  the  case  of  one  copy 
and  are  omitted  to  satisfy  page  limitations.  Now,  the  global 
optimization  procedure  is  easy  to  conceive,  as  follows: 

1 .  Begin  with  an  empty  list  of  local  optima. 

2.  Identify  and  list  one  phase  shift  vector  inside  each  sub- 
region  of  the  partition  of  the  c-dimensional  hypercube, 
the  feasible  region  of  phase  shift  vectors. 

3.  Pick  a  phase  shift  vector  from  the  list  in  Step  1.  Evaluate 
the  correspondingly  best  data  sequence  detection.  Re- 
estimate  the  phase  shift  vector  for  the  identified  data 
sequence  detection  and  re-evaluate  the  sequence 
detection  for  the  re-estimated  phase  shift. 

a.  If  the  two  sequence  detections  are  not  identical, 
the  sub-region  does  not  contain  a  local  optimum. 
Discard  the  results  over  this  sub-region. 

b.  If  the  two  sequence  detections  are  identical,  we 
have  a  local  optimum.  Insert  it  into  the  list  of 
local  optima. 

4.  Repeat  Step  3  for  each  phase  shift  vector  in  the  list 
created  in  Step  2. 
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5.  Search  the  list  of  all  local  optima  and  select  the  global 
optimum.  Output  the  corresponding  sequence  joint  ML 
detection  and  phase  shift  vector  estimate. 


The  difficulty  with  this  ML  solution  is  that  the  boundaries 
separating  the  sub-regions  are  curved  hyper-surfaces  as 
mentioned  above  and  are  computationally  cumbersome  to 
evaluate.  Therefore,  identifying  one  phase  shift  vector  inside 
every  sub-region  in  Step  2  is  an  impediment  to  implementing  the 
ML  solution.  To  keep  the  computations  at  an  acceptable  level, 
we  propose  the  following  algorithm  as  an  approximation.  Later 
on,  we  point  out  through  simulation  results  that  this  algorithm  is 
a  very  nearly  globally  optimal  for  the  application  domain  under 
consideration. 

IV.  THE  SIGNAL  COMBINING  ALGORITHM 
A  local  initial  point  for  the  phase  shift  parameter  vector  is 
obtained  by  maximizing  the  sum  of  components  in  the  objective 
function  corresponding  to  the  reference  signals  and  one 
particular  data  signal.  That  is,  let  (pij  represent  the  initial  phase 
estimate  for  the  y-the  copy  based  on  data  signal  i  only.  This  is 
determined  by  maximizing  as  follows: 

hi  =  max  y  \aij  cos{(pij  -  bij  +  6^) 

+  ^5yC0s((Piy-tiy)  .  (13) 

1=1  ' 


The  above  maximization  is  easily  carried  out  by  manipulating 
the  argument  as  a  linear  combination  of  cos^cpij)  and  sin((pjy) 
as  follows.  Let 


ijk  '  ijk’ 


and  ipk,  uniquely  determined  by 
cos{xpijk)  =7^,  and 

^ijk 

sin(i/)4fc)  = 


Then, 


-‘ijk 


^ik  ~  ^  ’  ^ijk^^^i.^ijk  ^ijk^f 
7=1 

c 

hi  =  max  >  QyfcCos((pyfc 

(pij,k 


(18) 

(19) 

(20) 

(21) 

(22) 


Let  k  =  k*  maximize  Then  the  y'-th  component  of  the  local 
reference  phase  vector  estimate  corresponding  to  data  signal  i  is 
given  by 

(Pij  =  ^liji,..  (23) 

The  initial  vector  of  phase  shifts  corresponding  to  data  signal  i  is 
given  by 

^i  =  (cpii.-.cpic),  i  =  (24) 

These  0;,  i  =  1, ■••,n  now  constitute  the  list  of  all  the  initial 
phase  vectors  replacing  the  list  in  Step  2  of  the  globally  optimum 
ML  signal  combining  solution  developed  in  Section  3.  This 
completes  the  development  of  our  algorithm  for  signal 
combining. 


1=1 


y=i 


=  ^  cos{(pijk)  cos(by  -  ^fc)  +  ^ N/Cos(t,y) 

+  sin((pyfc)  ^aij  sm{bij  -  6,,) 

+  '^Sijsm{tij)]. 


1=1 


Let 


V.  SIMULATION  EXPERIMENTS  AND  RESULTS 
We  conducted  extensive  simulation  experiments  and  compared 
the  performances  of  a  few  methods  for  phase  estimation, 
detection,  and  signal  combining.  The  first  of  these  is  the 
(14)  traditional  method  that  uses  only  the  reference  symbols  for  phase 
correction  and  evaluates  soft  bit  decisions  separately  in  each 
copy  and  then  adds  them  for  signal  combining.  The  second 
method  uses  the  ML  algorithm  for  phase  estimation  and  soft 
decision  evaluation  independently  for  each  copy  and  then  adds 
them  for  signal  combining.  Finally,  our  new  joint  ML-based 
signal  combining  algorithm  is  simulated.  The  above  algorithms 
are  applied  to  three  configurations  of  a  hop.  These  are  the  5-1-9 
mode  which  has  5  reference  symbols  and  9  data  symbols  per 
hop,  the  2-1-12  mode,  and  the  3-1-54  mode.  Recently,  we  obtained 
(1^)  virtually  identical  numerical  results  (as  from  the  above  joint  ML- 
based  algorithm)  with  a  more  efficient  algorithm  that  matches 
phase  shifts  before  joint  phase  estimation  and  detection. 
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BER  at  the  input  to  decoder  BER  at  the  input  to  decoder  BER  at  the  input  to  decoder 


Figure  1 .  Comparison  of  BERs  for  4  copies,  QPSK,  5+9 


Figure  4.  Comparison  of  BERs  for  4  copies,  8-PSK,  2+12 
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Figure  2.  Comparison  of  BERs  for  4  copies,  8-PSK,  5+9 
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Figure  5.  Comparison  of  BERs  for  4  copies,  QPSK,  3+54 
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Figure  3.  Comparison  of  BERs  for  4  copies,  QPSK,  2+12 


Figure  6.  Comparison  of  BERs  for  4  copies,  8-PSK,  3+54 
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We  experimented  with  quadrature  phase  shift  keying  (QPSK) 
and  8-PSK  (phase  shift  keying  with  eight  symbols  in  the 
eonstellation)  modulation  sehemes  eorresponding  to  2  and  3  bits 
per  symbol.  These  experiments  are  first  eondueted  with  an 
aeeurate  earrier  frequeney  and  a  uniformly  distributed  random 
phase  shift.  That  is,  in  this  ease,  there  is  no  phase  drift  over  the 
sueeessive  symbols  in  a  given  hop  and  the  loeation  of  referenee 
symbols  is  irrelevant.  In  another  set  of  experiments,  we  added  a 
phase  drift  for  eaeh  hop  to  simulate  a  frequeney  error.  In 
different  hops,  the  total  phase  drift  over  the  entire  hop  is 

uniformly  distributed  between  a  nominally  aeeeptable 

range  in  military  satellite  eommunieation  (MILSATCOM).  The 
referenee  symbols  also  experienee  phase  drift  and  so  their 
loeations  affeet  the  deteetion  perfomianee.  To  aeeount  for  this, 
the  loeations  of  referenee  symbols  are  randomized  over  the 
entire  hop.  Indeed,  their  loeations  over  multiple  eopies  of  the 
same  hop  are  independently  ehosen  at  the  transmitter  (but  they 
are  known  at  the  reeeiver). 

In  our  simulation  experiments,  all  the  data  symbols  in  every 
original  hop  are  generated  statistieally  independently  and  with 
equal  probabilities.  These  are  eonsidered  to  be  segments  of  one- 
third  rate  FEC  eneoded  and  interleaved  eode  words.  The  soft 
deeision  outputs  from  the  signal-eombined  deteetion  proeedure 
would  be  the  input  to  the  FEC  deeoder.  The  phase  shift 
experieneed  by  every  hop  is  simulated  as  independent  random 
numbers  in  the  range  [0,  Ztt).  All  the  eopies  of  a  hop  are 
subjected  to  AWGN  with  the  same  noise  power  spectral  density. 
In  order  to  illustrate  all  the  capabilities  of  our  joint  ML-based 
signal  combining  algorithm,  we  plot  results  for  the  experiments 
with  phase  drifts  and  randomized  locations  of  reference  symbols. 
Results  of  the  other  signal  combing  methods  in  the  plots  are 
based  on  no  phase  drift.  The  locations  of  reference  symbols  are 
irrelevant  when  there  is  no  phase  drift.  Results  for  combining 
four  copies  are  plotted  in  Figures  1  through  6.  We  conducted 
experiments  with  two  copies  also  and  some  summary  results  of 
these  experiments  are  given  in  Table  1 . 

The  horizontal  axis  represents  the  ratio  of  the  energy  in  each 
data  symbol  of  each  copy  to  the  noise  power  spectral  density. 
This  is  denoted  by  The  vertical  axis  represents  the  bit  error 

rate  resulting  from  the  different  signal  combining  algorithms. 
Soft  bit  decisions  from  the  signal-combined  detection  will  be 
input  to  the  FEC  decoder.  The  bit  error  rate  (BER)  required  at 
the  input  to  the  decoder  should  be  0.166  or  less  in  order  to  result 
in  a  BER  of  10“^  or  better  at  the  output  of  the  on-third  rate 
decoder.  The  range  of  —  over  which  the  results  plotted  is  the 

same  for  all  the  experiments  and  is  chosen  to  highlight  the  plots 
around  the  BER  value  of  0.166  for  each  plot.  Table  1  below  lists 
the  savings  in  dB  that  our  signal  combining  algorithm  offers 
over  the  power  that  would  be  required  for  the  traditional  signal 
combining  method  that  uses  only  reference  symbols  for  phase 
shift  correction,  for  the  case  of  BER  =  0.166.  The  columns  listed 
as  “Indep.  ML”  shows  the  power  savings  over  the  traditional 
method  for  the  case  of  using  the  ML  method  separately  for  each 


copy  before  adding  soft  decisions.  The  best  entry  in  the  table  is  2 
dB  for  four  copies  in  the  QPSK  2-1-12  mode. 


Table  1.  Power  savings  in  dB  over  the  traditional  method  at  BER 


=  0.166 


Mode 

Two  copies 

Four  copies 

Indep.  ML 

Joint  ML 

Indep.  ML 

Joint  ML 

QPSK  54-9 

0.0 

0.4 

0.7 

1.1 

8-PSK  5+9 

0.1 

0.1 

0.7 

0.7 

QPSK  2+12 

0.1 

0.9 

1.2 

2.0 

8-PSK  2+12 

0.1 

0.4 

1.3 

1.6 

QPSK  3+54 

0.0 

0.7 

1.0 

1.8 

8-PSK  3+54 

-0.1 

0.3 

0.9 

1.3 

VI.  CONCLUSION 

Low  rate  repeat  codes  are  useful  for  satellite  communication  in 
severely  contested  environments.  Frequency  hopping  is  often 
employed  in  these  systems  for  protection.  The  hypothetically 
best  gain  in  —  that  can  be  realized  by  retransmitting  once  after 

Nq 

an  original  hop  is  transmitted  is  3  dB  and  this  requires  coherent 
detection  which  is  not  feasible  in  FFI  systems.  In  our  simulation 
experiments,  we  observed  a  gain  of  only  2.3  dB  instead  of  the  3 
dB  for  the  QPSK  5-1-9  mode.  For  four  copies  in  the  same  mode, 
the  gain  is  only  4  dB  instead  of  6  dB.  In  severely  contested 
environments,  recovering  even  a  fraction  of  a  dB  with  a  better 
signal  combining  algorithm  is  very  valuable.  In  response  to  this 
requirement,  we  have  developed  a  new  signal  combining 
algorithm  that  is  based  on  the  maximum  likelihood  approach. 
Our  results  show  an  improvement  of  up  to  2  dB  over  the 
traditional  method  of  signal  combining  for  different  burst  modes. 
The  signal  combining  algorithm  is  very  robust  against  allowable 
local  oscillator  frequency  error  and  randomized  locations  for 
reference  symbols.  The  latter  feature  is  very  useful  in  preventing 
a  prospective  specialized  jammer  from  concentrating  its  power 
on  reference  symbols  only. 
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